Electromagnons and instabilities in magnetoelectric materials with non-collinear spin 
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We show that strong electromagnon peaks can be found in absorption spectra of non-collinear 
magnets exhibiting a linear magnetoelectric effect. The frequencies of these peaks coincide with the 
frequencies of antiferromagnetic resonances and the ratio of the spectral weights of the electromagnon 
and antiferromagnetic resonance is related to the ratio of the static magnetoelectric constant and 
magnetic susceptibility. Using a Kagome lattice antiferromagnet as an example, we show that 
frustration of spin ordering gives rise to magnetoelastic instabilities at strong spin-lattice coupling, 
which transform a non-collinear magnetoelectric spin state into a coUinear multiferroic state with 
a spontaneous electric polarization and magnetization. The Kagome lattice antiferromagnet also 
shows a ferroelectric incommensurate-spiral phase, where polarization is induced by the exchange 
striction mechanism. 

PACS numbers: 75.80.-|-q,71.45.Gm,76.50.-t-g,75.10.Hk 
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I. INTRODUCTION 

The recent renewal of interest in multiferroic materials 
led to discovery of many novel compounds where electric 
polarization is induced by ordered magnetic states with 
broken inversion symmetryiii^i^ The electric polarization 
in multiferroics is very susceptible to changes in spin or- 
dering produced by an applied magnetic field, which gives 
rise to dramatic effects such as the magnetically-induced 
polarization flops and colossal magnetocapacitancei^i^ 
Magnetoelectric interactions also couple spin waves to 
polar phonon modes and make possible to excite magnons 
by an oscillating electric field of light, which gives rise to 
the so-called electromagnon peaks in photoabsorptioni^ 

Electromagnons were recently observed in two 
groups of multiferroic orthorombic manganites, i?Mn03 
{R = Gd,Tb,Dy,Eui_:rY:r) and i?Mn205 (R = 
Y,Tb).^'^'-'^°-^^ Ferroelectricity in RMnOs appears in a 
non-collinear antiferromagnetic state with the cycloidal 
spiral ordering and the magnetoelectric coupling origi- 
nates from the so-called inverse Dzyaloshinskii-Moriya 
mechanism.— ii^iii'i^'i^ In Ref. it was noted that 
the same mechanism can couple magnons to photons 
and that an oscillating electric field of light can excite 
rotations of the spiral plane. However, the selection 
rule for the electromagnon polarization resulting from 
this coupling does not agree with recent experimental 
dat a-'^^d^d^i^Q and, moreover, the inverse Dzyaloshinskii- 
Moriya mechanism of relativistic nature is too weak to ex- 
plain the strength of the electromagnon peaks in i^MnOa. 

These peaks seem to originate from the exchange stric- 
tion, i.e. ionic shifts induced by changes in the Heisen- 
berg exchange energy when spins order or oscillatCfSi 
This mechanism explains the experimentally observed 
polarization of electromagnons. Since the Heisenberg ex- 
change interaction is stronger than the Dzyaloshinskii- 
Moriya interaction, it can induce larger electric dipoles. 
In Ref. [2l| it was shown that the magnitude of the spec- 



tral weight of the giant electromagnon peak in the spiral 
state of rare earth manganites is in good agreement with 
the large spontaneous polarization in the E-type antifer- 
romagnetic state^^ which has not been reliably measured 
yet but is expected to exceed the polarization in the spi- 
ral state by 1-2 orders of magnitude i^i^i 

From the fact that the mechanism that couples 
magnons to light in rare earth manganites is different 
from the coupling that induces the static polarization 
in these materials we can conclude that electromagnons 
can also be observed in non-multiferroic magnets. In this 
paper we focus on electromagnons in materials exhibit- 
ing a linear magnetoelectric effect, i.e. when an applied 
magnetic field, H, induces an electric polarization, P, 
proportional to the field, while an applied electric field, 
E, induces a magnetization, M. This unusual coupling 
takes place in antiferromagnets where both time reversal 
and inversion symmetries are spontaneously brokeui ^^i^^ 

It is natural to expect that when an electric field ap- 
plied to a magnetoelectric material oscillates, the induced 
magnetization will oscillate too. Such a dynamical mag- 
netoelectric response, however, requires presence of ex- 
citations that are coupled both to electric and magnetic 
fields. They appear when magnons, which can be excited 
by an oscillating magnetic field (antiferromagnetic reso- 
nances), mix with polar phonons, which are coupled to 
an electric field. Thus in materials showing a linear mag- 
netoelectric effect, for each electromagnon peak there is 
an antiferromagnetic resonance with the same frequency. 

This reasoning does not apply to all magnetoelectrics 
and the dc magnetoelectric effect is not necessarily re- 
lated to hybrid spin-lattice excitations. As will be dis- 
cussed below, in materials with collinear spin orders elec- 
tromagnons either do not exist or have a relatively low 
spectral weight. In this paper we argue that electro- 
magnons should be present in non-collinear antiferromag- 
nets showing strong static magnetoelectric response. As 
a simple example, we consider a Kagome lattice antifer- 
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FIG. 1: (Color online) The Kagome magnet with the KIT- 
Pite crystal structure, in which the ligand ions (open circles) 
mediating the superexchange between spins are positioned 
in a way that gives rise to a strong linear magnetoelectric 
response in the 120° spin state. Here, Ji and J2 denote, 
respectively, the nearest-neighbor and next-nearest-neighbor 
exchange constants, the solid arrows denote spins, while the 
empty arrows denote the shifts of the ligand ions. 



romagnet with the 120° spin ordering, shown in Fig. [T] 
Such an ordering has a nonzero magnetic monopole mo- 
ment, which allows for a linear magnetoelectric effect 
with the magnetoelectric tensor ofy = aSij for electric 
and magnetic fields applied in the plane of the Kagome 
lattice;^ A relatively strong magnetoelectric response 
was recently predicted for Kagome magnets with the 
KITPite crystal structure, in which magnetic ions are 
located inside oxygen bipyramidsi^ In this structure the 
oxygen ions mediating the superexchange in basal planes 
are located outside the up-triangles forming the Kagome 
lattice and inside the down-triangles or vice versa (see 
Fig. [1]), in which case magnetoelectric responses of all 
triangles add giving rise to a large magnetoelectric con- 
stant. 

This paper is organized as follows. In Sec.|TT]we analyze 
the symmetry of magnon modes and the magnetoelectric 
coupling in the Kagome lattice magnet with the KITPite 
structure and show that the dc magnetoelectric effect 
in this system is related to presence of electromagnon 
modes. The common origin of the dc and ac magneto- 
electric responses implies existence of relations between 
static and dynamic properties of magnetoelectric mate- 
rials, derived in Sec. IIIII In Sec. IIVI we discuss softening 
of (electro)magnons and the resulting divergence of the 
coupled magnetoelectric response. In Sec. |V] we discuss 
the transition from a magnetoelectric to a multiferroic 
state at a strong spin-lattice coupling and plot the phase 
diagram of our model system. In section IVII we discuss 
the importance of non-coUinearity of spins for dynamic 



magnetoelectric response and possible electromagnons in 
known magnetoelectric materials. In Sec. IVIII we con- 
clude. 



II. SYMMETRY CONSIDERATIONS 

The coupling of magnetic excitations in the Kagome 
magnet to electric field, resulting from the Heisenberg ex- 
change striction or any other non-relativistic interaction, 
can be found using the method outlined in Ref . [U . To 
simplify notation, we consider a single up-triangle, which 
has the same point symmetry as the whole Kagome lat- 
tice with the 120° spin ordering shown in Fig. [TJ The 
form of the magnetoelectric coupling is constrained by 
the 3z and nix symmetry operations:^^ 

ffmc - -7 I ^[(82-83) -(Si- S3)] 

+ ^[(Si.S3) + (S2-S3)-2(Si.S2)]|.(l) 



We then replace S.; by {S)ni + 6Si, where the unit vec- 
tors (ni,n2,n3) = ^-^i - iy, - iy, describe 

the 120° spin ordering in the xy plane and 6Si _L is 
the oscillating part, which is the superposition of the or- 
thogonal magnon modes in the triangle (the zero wave 
vector magnons for the Kagome lattice) : 



<5S, = (^-V'™ + {S)P 



(2) 



where a — 0,x,y labels the magnon, 

c^o. = (1,1,1), 



(3) 



are the out-of-plane components of the magnons and 



rii are the in-plane components (see Fig. [2]), 
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The singie-magnon excitation by an electric field is 
described by the terms linear in dS, while the terms 
quadratic in (5S give rise to the photoexcitation of a 
two-magnon continuum. Since spins order in plane, the 
polarization oscillations are induced by the in-plane os- 
cillations of spins and the coupling of electric field to 
magnons, obtained from Eq.((T]), has the form: 



-9e {qxEx + qyEy) , 



(5) 



where gs = §7(5'). This magnetoelectric coupling is 
only nonzero in the magnetically ordered state with bro- 
ken time reversal symmetry, which is why the coupling 
constant is proportional to (S). 
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FIG. 2: (Color online) The magnon modes in a triangle with 
the 120° spin ordering. The thin arrows indicate the direc- 
tions of average spins, while the thick arrows show the in- 
plane components of the magnon modes. 



We note that Eq. ([5]) can also be obtained by a conven- 
tional symmetry analysis of magnetic modes. The vector 
tpQ and, hence, the corresponding amplitude, go, forms 
a one-dimensional representation Fi, while il'y) ^'Hd 

qx 



Qy 



form a two-dimensional representation F3 (see Ta- 



ble IJ). The direct product, F2 x F3, where F2 is the 
symmetry of the spin ordering, transforms as the dou- 
blet of the in-plane components of the electric fields, 
which leads to Eq.([5]). This general symmetry analysis is, 
however, insensitive to the microscopic mechanism of the 
magnetoelectric coupling, whereas the derivation staring 
from Eq.ijT]) only applies to non-relativistic mechanisms. 
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TABLE L The transformation properties of several irreducible 
representations of the space group of the Kagome lattice and 
time reversal operation T. 

Table |T] shows that the coupling to the in-planc com- 
ponents of magnetic field has the form, 



9h {qxHx + qyHy 



(6) 



while the Hamiltonian describing magnon modes with 
zero wave vector in the Kagome layer is, 



H{p,q) 



vl 



1 
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+ 2 ('^'+'^')- 
(7) 

If we consider, for example, the microscopic spin Hamil- 
tonian describing the antiferromagnetic nearest-neighbor 
and next-nearest-neighbor Heisenberg exchange interac- 
tions with the exchange constants, respectively, 3\ and 



J2 and the easy plane magnetic anisotropy A, 

= Ji^S,.S, + J2 S.-S,-f |^^(5ff , (8) 



for which the 120° spin ordering shown in Fig. [T] is a 
classical ground state,— we get 



TO-i = [6(J + J') + A](5)2, m-i = A(5)2, 
Ko = 0, K = 3 ( Ji + J2) . 



(9) 



The linearized equations of motion for spins in applied 
electric and magnetic fields are obtained from Eqs.®,®, 
and ([7]), if we impose the commutation relations for the 
amplitudes of the in-plane and out-of-plane parts of (SS: 



(10) 



From these equations we find the frequencies of the three 
magnon modes with zero wave vector: cjg — K^m^^ — 
and iJl=ij?y = Km-^ =3{Ji + J2) A{Sf. 

Minimizing the spin energy with respect to q^ and qy 
in the static limit, we obtain an effective magnetoelectric 
coupling, 



HaiC — —a [Hr^Ej. + HyEy) . 



(11) 

where the magnetoelectric coefficient a = SEgn ^ Pu^- 
thermore, the ga;-mode can be excited by both electric 
and magnetic field oscillating with the frequency tOx in 
the direction parallel to the x axis, while the g^-mode 
can be excited by both Ey and Hy, which shows that the 
static linear magnetoelectric effect in this non-coUinear 
magnet is related to the presence of elcctromagnon and 
antiferromagnetic resonance peaks with equal frequencies 
in the optical absorption spectrum. 



III. RELATIONS BETWEEN STATIC AND 
DYNAMIC MAGNETOELECTRIC RESPONSE 

The common origin of the static and dynamic mag- 
netoelectric response of non-coUinear magnets leads to 
quantitative relations between dc susceptibilities and 
spectral weights of peaks in the optical absorption spec- 
trum. These relations simplify when the coupling of 
magnons to electric field is mediated by a single optical 
phonon. The description of magnons in terms of con- 
jugated coordinates and momenta is very convenient for 
derivation of these relations, since the coupled magnon 
and optical phonon are in this approach just a pair of 
coupled oscillators: 
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(12) 
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where {Qxi Px) and {Qy, Py) are the coordinates and mo- 
menta of the optical phonons coupled to magnons, which 
also form a two-dimensional representation. 

The magnetoelectric response of such a system is easy 
to calculate. The result can be expressed in terms 
of observable quantities, such as the 'dressed' magnon 
and phonon frequencies, Wmag and tjph, and the spec- 
tral weights of the magnon and phonon peaks excited by 
an electric and magnetic field. We denote the spectral 
weight of the electromagnon peak by 



5,1, = 8 ldu;cjx'^{cj), 



(13) 



where Xc i^) the imaginary part of the dielectric ac 
susceptibility, while the spectral weight of the antiferro- 
magnetic resonance is. 



(14) 



where Xmi^) ^he imaginary part of the magnetic ac 
susceptibility. The integration in Eqs. p^ and Eq. p^ 



goes over an interval of frequencies around Wmag- The two 
spectral weights for the phonon S'^, and are defined 
in a similar way. We assume that the magnon and phonon 
peaks are sufficiently narrow and can be separated from 
each other. The four spectral weights satisfy a relation, 

(15) 

following from the fact that an electric field only interacts 
with the 'bare' phonon, while a magnetic field is only 
coupled to the 'bare' magnon. 

The relations between the dc and ac magnetoelectric 
responses of the coupled spin-lattice system have the 
form. 



*-^mag^mag 
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(16) 
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where Ae(A/i) is the increase of the real part of the 
dielectric constant (magnetic permeability) at zero fre- 
quency resulting from the magnon and phonon peaks 
(we use the Gaussian units). The first two equations 
are the Kramers-Kronig relations for the real and imag- 
inary parts of dielectric and magnetic susceptibilities, 
while the last relation follows from equations of motion. 
Combining Eqs.(fT5| and (fT6|) we can express the ratio of 
the spectral weights of the electromagnon and the anti- 
ferromagnetic resonance through the ratio of the static 
magnetoelectric constant a and magnetic susceptibility 

Xm = Xm(^ = 0): 



5. 



cH 
^mag 



1 



ph 



(17) 



For w,lg <C Wph, the ratio of the spectral weights is just 
the square of the ratio of the dc magnetoelectric constant 
and magnetic susceptibility. 

Due to the spin-lattice coupling, the 'dressed' magnon 
frequency, ciJniag, is lower than its 'bare' value, (as- 
suming that the 'bare' magnon frequency is smaller than 
the 'bare' phonon frequency). As the spin-lattice cou- 
pling increases, Wmag vanishes at a critical value of the 
coupling. According to Eg. p^ . this results in the simul- 
taneous divergency of e, fi and a, indicating an insta- 
bility towards a multiferroic phase, which is both ferro- 
electric and ferromagnetic. Another manifestation of this 
instability is the fact that as the spin-lattice coupling ap- 
proaches the critical value, the magnetoelectric constant 
a tends to its upper bound equal the geometric mean 
of the dielectric and magnetic susceptibilities, y^XoXm, 
imposed by the requirement of stability with respect to 
applied electric and magnetic fieldsi^ 



IV. MAGNON SOFTENING 

To study the transition from the magnetoelectric state 
of the Kagome magnet to the multiferroic state in more 
detail, we consider a simple microscopic model, in which 
positions of magnetic ions are fixed, while ligand ions 
are allowed to move. The spin-lattice coupling originates 
from the dependence of the exchange constants on dis- 
placements of ligand ions mediating the superexchange in 
the direction perpendicular to the straight line connect- 
ing two neighboring spins. We denote the positions of the 
three ligand ions outside up-triangles by Ui, U2, and U3, 
while the position of a ligand ion inside a down-triangle 
is denoted by v. Then, for example, the exchange con- 
stant for the spins Si and S2 is Ji -I- J( (^3)^^ , while for 
the spins S4 and S5 it is Ji + J[vy (see Fig. [T]). Further- 
more, we assume that phonons are dispersionless and the 
lattice energy for a pair of the up- and down-triangles is. 



Ulat - 




(18) 



where K is the spring constant. 

For the 120° structure shown in Fig. [3]^ a), the magne- 
toelectric constant a in Eq. lTTI) is given by 



{1- g)2{Ji + J2)Kv' 



(19) 



where Q = — 2e is the charge of the oxygen ion, (/i) — 
2^.b{S) is the average magnetic moment on each site, v 
is the unit cell volume, and 
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15 {J[{S)f 



8 (Ji + J2)K 
is the dimensionless spin-lattice coupling constant. 



(20) 
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An estimate for the magnetoelectric constant, a ~ 
10~^, for the model parameters appropriate for the KIT- 
Pite structure (S = 2, Ji - 3meV, K - 6eV • A"^ 
i - 3.5A" , and 1; = I77A ) agrees well with the result 
of ab initio calculationsi^ Furthermore, Xm ^ 2 • 10^^, 
so that (^^^ 25. Thus, the electromagnon peak in 

KITPite should be much stronger than the antiferromag- 
netic resonance peak, which is also the case for rare earth 
manganites with a spiral ordering. 

At 5 = 1, the magnetoelectric constant diverges and 
so do the dielectric and magnetic susceptibilities: 



1 



1 



(21) 



Since in our model there are two polar phonons coupled 
to a magnon with a given polarization (one in the up- 
triangle and another in the down-triangle) , the magneto- 
electric constant comes close to its upper bound but does 
not reach it at g = 1: 



0.985. 



(22) 



3=1 



Surprisingly, the softening of the magnon mode, 
which is not coupled to polar lattice distortions, occurs 
at a lower value < \. Since loq — the softening 
in this case means vanishing velocity of the go-mode. 
The velocity vanishes, because away from the F-point 
in the magnetic Brillouin zone magnons with different 
symmetry become mixed and the go-mode is coupled to 
the electromagnon modes. As the spin-lattice coupling 
grows and the electromagnon frequency decreases, the 
lowest-frequency magnon branch is pushed down, which 
ultimately reduces the velocity of the go-mode to zero. 



V. MAGNETOELASTIC INSTABILITIES 

Although KITPite, for which g ~ 0.05, is far away from 
the instabilities discussed in the previous section, it is in- 
teresting to study behavior of magnetoelectric materials 
when the spin-lattice coupling becomes strong, in partic- 
ular, in view of the dramatic magnetoelectric effects re- 
cently observed in multiferroics. As the magnetoelectric 
constant becomes large close to the transition between 
magnetoelectric and multiferroic states, it is important 
to understand possible scenarios of such a transition. 

In this section we present the analytical and numeri- 
cal study of the phase diagram of the KITPite layer for 
strong spin-lattice couplings. In particular, we show that 
none of the continuous transitions involving magnon soft- 
ening, discussed in Sec. IIVI actually takes place, as the 
strong spin-lattice coupling makes the frustrated Kagome 
magnet unstable towards a first-order magnetoelastic 
transition that relieves the frustration. This frustration- 
driven instability is similar to the one found in spinels, 
where a collinear ordering of spins appears together with 



(a) 



(b) 



(c) 




FIG. 3: (Color online) The minimal-energy spin configura- 
tions of the Kagome magnet for three different values of the 
spin-lattice coupling: (a) the 120°-state with zero wave vec- 
tor, (b) the incommensurate ferroelectric state, and (c) the 
collinear multiferroic state. The short solid arrows show the 
spin directions, while the short empty arrows denote the shifts 
of the ligand ions in these states. The long solid and empty 
arrows show the direction of, respectively, the spontaneous 
magnetization and polarization. 
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a lattice deformationi^i^i^i^ We show that the trans- 
formation of a non-collinear magnetoelectric state into 
a colhnear multiferroic state can involve two transitions 
and an intermediate phase, which is ferroelectric but not 
ferromagnetic. 

To understand the origin of magnetoelastic instabili- 
ties at strong spin-lattice coupling, we first consider a 
single up-triangle and integrate out the lattice degrees of 
freedom, (i = 1,2,3). Then the total energy of the 
triangle takes the form of an effective spin Hamiltonian 
with quadratic and bi-quadratic interactions: 



A 



E 



JiSi ■ So 



{J'lf 

2K 



(23) 



The bi-quadratic interactions favor coUinear spins and 

fo": 9 = = ^^^5 > i the lowest-energy 

spin configuration is a collinear state of the ffj type 
(the spins lie in the lattice plane). Similarly, an effec- 
tive spin Hamiltonian for a down-triangle, where the ex- 
change along all bonds is mediated by a single ligand ion 
located inside the triangle, has the form. 



E 



JiSi ■ Sj „ (Si • Sj) 



2K 



+^ E (s.-s,)(s,.s,). 



(24) 



In this case the critical coupling is lower: g — ||. 
Due to the inequivalence of the up- and down-triangles 
in the KITPite structure the transition from the 120°- 
state shown in Fig. [3^ to a fully collinear state shown in 
Fig. [3t goes in two steps via an intermediate state, where 
only the spins in the down-triangles are collinear while 
the spins in the up-triangles are still non-collinear (see 
Fig.Eb). 

The transition to the intermediate state does not fully 
lift the frustration: the number of combinations of the 
down-triangles with collinear spins and up-triangles with 
the 120°-angle between spins grows exponentially with 
the system size. In our model this degeneracy is re- 
moved by the next-nearest-neighbor interactions between 
spins, which select the state shown in Fig. [Sb. The next- 
nearest-neighbor interactions in the vertical direction in- 
duce a small spin canting, as a result of which spins in 
down-triangles are not strictly collinear and the angles 
between spins in up-triangles deviate from 120° by the 
angle ip Furthermore, the nearest-neighbor in- 

teractions in the remaining two directions are frustrated 
in the commensurate spin state with the wave vector 

Q = 2fc(l,0) + k (^i, where k ^ ^ {a being the 

lattice constant). This frustration is lifted, when the 
spin ordering becomes incommensurate with the lattice: 
k^ ^+6, where 6 (x ^. 

3a ' Ji 

The incommensurate state has zero net magnetization, 
but its electric polarization is nonzero. For the state 
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FIG. 4: (Color online) Plotted is the zero-temperature phase 
diagram of the Kagome layer for large values of the spin- 
is OHf))! _ For relatively 
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lattice coupling g 
small ^ the magnetoelectric (ME) phase with the 120° spin 
ordering [see Fig.[3ja)] undergoes a first-order transition into 
the incommensurate (IC) spin state, which is ferroelectric 
(FE) [see Fig.[3|b)], as the coupling constant g increases. Fur- 
ther increase of g results in the transition into the collinear 
multiferroic (MF) phase [see Fig. EJc)] , with the parallel spon- 
taneous polarization and magnetization, P || M. For larger 
■j^, the ME state undergoes a direct transition into the fully 
collinear MF state. Also plotted are the dot-dash line, where 
the go-mode in the ME state softens and the dotted line, at 
which the magnetoelectric response of the ME state diverges 
(both phenomena do not occur because of the intervening 
first-order transitions to the FE and MF states). 



shown in Fig. [H^b) the polarization vector is parallel to 
the y axis. This polarization originates not from the in- 
verse Dzyaloshinskii-Moriya mechanism, which makes in- 
commensurate spiral states in e.g. i?Mn03 ferroelectric, 
but from the fact that all bonds connecting parallel spins 
are parallel to each other, which forces the ligand ions in 
all down-triangles to shift in the same direction, since 



Vx = (S3 • S5 - S3 • S4) , 



(25) 



rj^ (2 (S4 • S5) - S3 • S4 - S3 • S5 



AT /' 

where r] — \/ and the labeling of spins is the same as 

in Fig. [TJ For the state shown in Fig. El^b), where bonds 
connecting (nearly) parallel spins are oriented along the 
X axis, Vx = 0, while Vy and hence the polarization Py is 
nonzero. 

We note that the coupling of exchange interactions 
to strains, which gives rise to magnetoelastic transitions 
in frustrated spinels f^i^ii^ also favors the ferroelectric 
state. The absence of inversion symmetry in the KIT- 
Pite layer allows for the piezoelectric coupling, 



(26) 
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where Uy is the strain tensor, so that the hgand displace- 
ment, V, the electric polarization, P, and the strains are 
coupled to each other. 

The zero-temperature phase diagram of the Kagome 
layer with g and ^ along the horizontal and vertical 

axes is shown in Fig. [4] For small the incommensu- 
rate ferroelectric (IC FE) state, discussed above, inter- 
venes between the magnetoelectric (ME) state with the 
120° spin ordering [see Fig. ^a,)] and the fully coUinear 
multiferroic (MF) state shown in Fig. [3]Jc), in which the 
spontaneous polarization and magnetization are parallel 
to each other, P || M. As the ratio ^ grows, the interval 
of the coupling constant g where the intermediate state 
is stabilized shrinks and above the tricritical point the 
ME state undergoes a direct transition into the coUinear 
MF state along the critical line ^ = |g — 1. Also plotted 
are the dot-dash line, at which the go-mode would soften 
and the dotted line, at which a, Xej and Xm would di- 
verge, if the ME state would survive at strong spin-lattice 
couplings. 

VI. DISCUSSION 

We showed that the static magnetoelectric response of 
non-coUinear antiferromagnets can be related to hybrid 
magnon-phonon modes coupled to both electric and mag- 
netic fields. Such magnetoelectric materials are analogs 
of displacive ferroelectrics the dielectric response of which 
is governed by optical phonon modes. If spins in an 
ordered state are coUinear, the exchange striction can- 
not couple an electric field to a single magnon, as the 
expansion of scalar products of parallel or antiparallel 
spins begins with terms of second order in 6S, which 
give rise to photoexcitation of a two-magnon continuum 
(the so-called "charged magnons"— ) . Electromagnons in 
coUinear magnets can still originate from mechanisms in- 
volving relativistic effects, such as the exchange striction 
induced by the antisymmetric Dzyaloshinskii-Moriya in- 
teraction, which is proportional to the vector product 
of two spins. In 3c? transition metal compounds such 
couplings are weak compared to the exchange striction 
driven by the Heisenberg exchange, so that the spectral 
weight of electromagnons in coUinear magnets should be 
relatively low. 

Magnetoelectric materials with coUinear spin orders 
may rather be analogs of 'order-disorder' ferroelectrics 
with the static magnetoelectric response originating from 
thermal spin fiuctuations. Cr203 seems to be an example 
of such a material: its magnetoelectric coefficient passes 
through a maximum below Neel temperature and then 
strongly decreases when temperature goes to zero and 
spin fluctuations become suppressed i^^*^ 

We note that the rotationally invariant coupling Eq. ^ 
may also originate from purely electronic mechanisms, 
such as the polarization of electronic orbitals induced 
by a magnetic orderingi^^iHiSliSSiiS Ab initio calcula- 
tions suggest that in rare earth manganites the electronic 



mechanisms of magnetoelectric coupling are as important 
as the exchange striction.— On the other hand, the in- 
crease of the spectral weight of the electromagnon peaks 
in iJMnOs below the spiral ordering temperature oc- 
curs largely at the expense of the strength of the optical 
phonon peak at ~ lOOcm^^, suggesting the dominant 
role of the spin-lattice coupling. ^^'^^ If electronic mech- 
anisms dominate and an electromagnon gets its spectral 
weight from frequencies much higher than those of opti- 
cal phonons, Eq. (fT6l) should to be modified in an obvious 
way, while Eq. (fT7|) . where should be replaced by 0, 
is still valid. 

We note that the non-coUinearity of spins by itself does 
not guarantee strong magnetoelectric effect and electro- 
magnon peaks - the crystal structure is equally impor- 
tant. Thus, in the layered Kagome antiferromagnet, the 
iron jarosite KFe3(OH)6(S04)2,^^ which has the spin or- 
dering shown in Fig. [1] the ligand ions are located out- 
side of both up- and down-triangles, which cancels the 
magnetoelectric effect due to the Heisenberg exchange 
striction. The cancellation also occurs in triangular mag- 
nets with the 120° spin ordering, as they contain three 
different spin triangles, such that spins in one triangle 
are rotated by ±120° with respect to spins in two other 
trianglesSS, (more generally, the linear magnetoelectric ef- 
fect can only be induced by a spin ordering with zero wave 
vector). We note, however, that the lattice trimerization 
in hexagonal manganites^ makes the three types of spin 
triangles inequivalent and destroys the cancellation. This 
can be also seen from the symmetry properties of the A1.2 
and Bi 2 phases of hexagonal manganites^ allowing for 
the magnetoelectric term E^Hy ~ EyH^ in the Ai-phase, 
which has a toroidal moment, and the term E^Hx+EyHy 
in the A2-phase, which has a magnetic monopole mo- 
ment. Whether electromagnons in these phases can be 
observed, depends on the magnitude of the trimerization 
and remains to be explored4^ 



VII. CONCLUSIONS 

In conclusion, we showed that magnets with non- 
coUinear spin orders resulting in a linear magnetoelec- 
tric effect may also show electromagnon peaks in optical 
absorption spectrum. While electromagnons should be 
present in many non-coUincar magnets, the specific fea- 
ture of magnetoelectric materials is that some magnon 
modes can be excited by both electric and magnetic 
fields, i.e. electromagnons are also antiferromagnetic res- 
onances. We derived a simple relation Eq. (|17p between 
the ratio of the spectral weights of the electromagnon 
and antiferromagnetic resonance peaks and the ratio of 
the static magnetoelectric constant and magnetic suscep- 
tibility, which can be used to estimate the strength of 
electromagnon peaks on the basis of dc measurements. 

To make our consideration more specific, we consid- 
ered a Kagome lattice magnet with the KITPite struc- 
ture, where the ligand ions are positioned in a way that 
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gives rise to a relatively strong linear magnetoelectric 
effect.^'* Using the symmetry analysis we identified the 
magnon modes that are coupled to both electric and mag- 
netic fields and give rise to the linear magnetoelectric ef- 
fect. We showed that the softening of these modes at a 
strong spin-lattice coupling results in the divergence of 
the magnetoelectric constant as well as of magnetic and 
dielectric susceptibilities, signaling the instability of the 
magnetoelectric state towards a multiferroic state with 
spontaneously generated P and M. However, the de- 
tailed study of the phase diagram of this model revealed 
that the electromagnon softening does not actually take 
place, since the first-order transition to the coUinear mul- 
tiferroic state occurs at a lower value of the spin-lattice 
coupling. In some region of model parameters a ferroelec- 
tric incommensurate-spiral phase intervenes between the 
magnetoelectric and multiferroic phases. While in known 
spiral magnets, ferroelectricity is likely induced by the 



inverse Dzyaloshinskii-Moriya mechanism, in our model 
it results from the stronger exchange striction mecha- 
nism due to the coUinearity of spins in half of the tri- 
angles. These magnetoelastic instabilities are typical for 
frustrated magnets, where non-coUinear spin orders usu- 
ally occur. We also discussed the possibility to observe 
electromagnons in known magnetoelectric materials. We 
hope that our study will stimulate experimental work in 
this direction. 
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